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Abstract 

Low- and intermediate mass Higgs bosons decay preferably into fermion pairs. The 
one-loop electroweak corrections to the respective decay rates are dominated by a 
flavour-independent term of OiGpm^). We calculate the two- loop gluon correction 
to this term. It turns out that this correction screens the leading high-mt behaviour 
of the one-loop result by roughly 10%. We also present the two- loop QCD correc- 
tion to the contribution induced by a pair of fourth-generation quarks with arbitrary 
masses. As expected, the inclusion of the QCD correction considerably reduces the 
renormalization-scheme dependence of the prediction. 

1 Introduction 

One of the great puzzles of elementary particle physics today is whether nature makes use of 
the Higgs mechanism of spontaneous symmetry breaking to generate the observed particle 
masses. The Higgs boson, H, is the missing link sought to verify this concept in the Standard 
Model. Many of the properties of the Higgs boson are fixed, e.g., its couplings to the gauge 
bosons, gvvH = 2 b ^G l p 2 M v (V = W, Z), and fermions, gjj H = 2 1 ^ 4 G 1 J 2 rrif, and the vacuum 

expectation value, v = 2 _1//4 G F 1 ^ 2 m 246 GeV. However, its mass, M#, and its self-couplings, 
which depend on Mh, are essentially unspecified. 

The failure of experiments at LEP 1 and SLC to detect the decay Z — > ffH has ruled 
out the mass range Mh < 63.8 GeV at the 95% confidence level [[L]]. At the other extreme, 
unitarity arguments in intermediate-boson scattering at high energies || and considerations 
concerning the range of validity of perturbation theory || establish an upper bound on Mh 

at ~ 1 TeV in a weakly interacting Standard Model. 

The Higgs-boson discovery potential of LEP 1 and SLC is almost exhausted @] . Prior to 
the advent of the LHC, the Higgs-boson search will be restricted to the lower mass range. 
With LEP 2 it should be possible to find a Higgs boson with Mh < 100 GeV when high 
energy and luminosity can be achieved ||. A possible 4-TeV upgrade of the Tevatron might 
cover the Mh range up to 120 GeV or so ||. At an e + e~ linear collider operating at 300 GeV, 
50 fb _1 luminosity and a 6-tagging efficiency of 50% would be sufficient to detect a Higgs 
boson with Mh < 150 GeV in the u, + u~bb channel 0. 
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Below the onset of the W + W~ threshold, the Standard- Model Higgs boson is relatively 
long-lived, with T H < 100 MeV, so that, to good approximation, its production and decay 
processes may be treated independently. The low-mass Higgs boson, with Mjj < Mz, decays 
with more than 99% probability into a fermion pair ||. With Mh increasing, the W + W~ 
mode, with at least one W boson being off shell, gradually gains importance. Its branching 
fraction surpasses that of the r + r~ mode at Mh ~ 115 GeV and that of the bb mode at 
Mh ~ 135 GeV |§. In the near future, however, Higgs-boson searches will rely mostly on 
the // modes. 

Quantum corrections to Higgs-boson phenomenology have received much attention in the 
literature; for a review, see Ref. ||. The experimental relevance of radiative corrections to 
the // branching fractions of the Higgs boson has been emphasized recently in the context 
of a study dedicated to LEP 2 ||. Techniques for the measurement of these branching 
fractions at a \fs = 500 GeV e + e~ linear collider have been elaborated in Ref. 
QCD corrections to the H — > qq decay rates are most significant numerically 
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In the 



approximation m g -C Mh, they are known to O (aj) |JL2|| . The theoretical uncertainty related 
to the lack of knowledge of the terms of O ^a 2 m 2 /M H ) and O (af ) is presumably small [ 13 
The bulk of the QCD corrections is attributed to the running of m q 



up to scale Mh- In the 
case of the bb mode, the QCD correction relative to the Born approximation implemented 
with the pole mass ranges between —53% and —63% for M H between 60 GeV and 2M W 0. 

The leading high-M# correction to the H —*■ ff decay widths is flavour independent. 
The one-loop term, of 0{GpM H \ was first obtained by Veltman |TJ|; it is positive and 

1 TeV. Recently, the two-loop O (GpM H ) term has been found [ |15[| ; 

400 GeV. 



reaches 11% at Mh 

it is negative and exceeds in magnitude the O (GpM H ) term already at Mh 
The leading contributions due to new heavy fermions are also independent of the final-state 
flavour; at one loop, they are positive and increase quadratically with the heavy-fermion 
masses [ I6|| . The full one-loop electroweak corrections to the H —> ff decay widths are 
now well established |17jl8||. They consist of an electromagnetic and a weak part, which 
are separately finite and gauge independent. The electromagnetic part emerges from the 
one-loop QCD correction jO]] by substituting aQ 2 ? for a s Cp, where Qf is the electric charge 
of / and C F = (N 2 - 1) /(4iV c ), with N c = 3. For M H < 2M W , the weak part is well 
approximated by [0 
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M^r/Mz, Kb = 1, and Kf = 7 for all other flavours, except for top. 



The tt mode will not be probed experimentally anytime soon and we shall not be concerned 
with it in the remainder of this paper. 

Throughout this paper, we adopt the so-called modified on-mass-shell (MOMS) scheme 
T9| , which emerges from the ordinary electroweak on- mass-shell scheme [0,0 by elimi- 



nating a in favour of Gp by virtue of the relation Gp = (na / y/2s 
Here, Ar embodies the non-photonic correction to the muon decay width. 



approximation of the MOMS scheme, one has |22 
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where N c = 1 for lepton flavours, and the weak correction is implemented by including the 
overall factor (1 + A weak ). 

Equation ([!]) has been obtained by putting Mh — mf — (/ ^ t) in the expression 
for the full one-loop weak correction. It provides a very good approximation for / = r up 
to Mh ~ 135 GeV and for / = b up to Mh ~ 70 GeV, the relative deviation from the 
full weak correction being less than 15% in each case. From Eq. ([]]) it is evident that the 
dominant effect is due to virtual top quarks. In the case / ^ b, the m t dependence is carried 
solely by the renormalizations of the wave function and the vacuum expectation value of 
the Higgs field and is thus flavour independent. These corrections are of the same nature as 



those considered in Ref. ||16|| . For f = b, there are additional m t dependent contributions 
from the bbH vertex correction and the 6-quark wave-function renormalization. Incidentally, 
they cancel almost completely the universal m t dependence. It is amusing to observe that 



a similar situation has been encountered in the context of the Z — > // decays p3| . In 



summary, the universal virtual-top-quark term will constitute the most important part of 
the weak one-loop corrections to Higgs-boson decays in the near future. In this paper, we 
shall present the two-loop gluon corrections to this term. 

This paper is organized as follows: In Sect. 2, we shall consider the universal contribution 
to the H —> f f decay rates induced by a pair of quarks with arbitrary masses and evaluate 
its QCD correction adopting the on-shell definition of quark mass. In Sect. 3, we shall 
repeat this calculation in the modified minimal-subtraction (MS) scheme [§4] and compare 
the result with the one of Sect. 2 in order to estimate the theoretical uncertainty related to 
the arbitrariness of the definition of quark mass. Section 4 contains the numerical analysis. 
Our conclusions are summarized in Sect. 5. 



2 Two-loop results 

In this section, we shall present the QCD correction to the shift in T (^H — > ff^j induced by a 
pair of quarks, (U, D), with arbitrary masses. For simplicity, we shall assume that U and D 
do not mix with /. Here, we shall adopt the on-shell definition of quark mass. As explained 
in Sect. 1, in the MOMS scheme, such corrections reside inside the renormalizations of the 
wave function and the vacuum expectation value of the Higgs field. The relevant part of 

^weak is 

i=-?=i-Mi„(M;), (3) 

where n^^fs) and Uhh(s) are the unrenormalized self-energies of the W and Higgs bosons, 
respectively, evaluated at four-momentum squared s. In the following, we shall write down 
only the (U, D) contributions to the quantities under consideration. 

For completeness, we shall first review the one-loop results. In dimensional regularization, 
they read 
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4 — 2e is the dimensionality of space-time, and £q = In (ttiq / fi 2 ^J , with fi being the 
't Hooft mass. Here and in the following, we suppress terms containing 7^ — ln(47r), where 
is Euler's constant. These terms may be retrieved by substituting /i 2 — > 4ire~ lE fi 2 . When 
the on-shell definition of mass is employed, fi will drop out in the expressions for physical 
quantities and so will these terms. In the MS scheme, these terms are subtracted along with 
the poles in e. Inserting Eqs. (f|||) in Eq. (Q), one obtains [0 
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In the limit mu^mr, 3> Mjj/2, this simplifies to (TB 
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Setting my = m t and mc = 0, one recovers the mt- dependent term of Eq. (|l|). Light quarks, 
with mu,rriD -C Mh/2, decouple from 5. 

We now proceed to two loops. The QCD correction to the Higgs-boson self-energy arises 
from the Feynman diagrams depicted in Fig. [I]. Apart from the two-loop gluon-exchange 
diagrams [see Fig. |l](a)], one needs to include the diagrams where the one-loop quark mass 
counterterm, 5mg, is inserted into the vertices and propagators of the one-loop seed diagram 
[see Fig. 0(b)]. 

In the on-shell scheme of mass renormalization, 5mq is adjusted in such a way that the 
pole of the renormalized propagator appears at the renormalized mass. Specifically, 
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is the quark self-energy due to the exchange of one virtual gluon. As before, the combination 
7# — ln(47r) has been absorbed into a redefinition of the 't Hooft mass. A straightforward 



calculation yields p6 
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Using a notation consistent with Ref. [26|, we find 
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where Ai and Yi are the divergent constants that have been introduced in connection with 
the gauge-boson vacuum polarizations at (D(aa s ) [pTjpSfl and Hi is a finite function. The ex- 
pressions for X\ and Y% depend on the regularization scheme. In dimensional regularization, 
they readP] 
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For r < 0, 



#i(r) 



6Li 3 (r 2 ) - 3Li 3 (ri) + 8/ (Li 2 (r 2 ) - Li 2 (rt)) + Af(-3f 
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+ ^(-C>,- + 2 + ^]+3/(-3/- + £ 



Li 2 (rVi - Li 2 (rM + /(-3/ + 2g + 4fc) 



^S + 3C(2)(r-3) + l(14) 
r 4 
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where Li 2 and Li 3 are the dilogarithm and trilogarithm f3(| , respectively, r± = y/1 — r ± 
y/—r, f = lnr + , g = ln(r + — r_), and h = ln(r + + r_). A table of handy transformation 
rules for the analytic continuation in r is available from Ref. []31 |. For r > 1, fZj develops 
an imaginary part, which, by Cutkosky's rule |3"2|| , is related to the cuts of the two-loop 
amplitudes. This provides the opportunity for a nontrivial check of the calculation, since 
this imaginary part is related to the well-known 0(a s ) correction to V (H — > qq). In fact, 
for Mu > 2m q , one verifies that 
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1 As a consequence of a misprint in Ref. P9| , the term of Y% involving £(2) occurs in Ref. fl28| with a wrong 
prefactor. Fortunately, this is inconsequential for the physical results of Ref. p8| because the Y\ terms cancel 
exactly among themselves. 
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where £qed is given by Eq. (3.9) of Ref. |p~8[ . 



The gluon correction to the (U, D) contribution to il W vi/(0) may be found in Ref. [[33 
where the on-shell definition of quark mass is employed. The result may be written as 
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Note that F(u,d) = F(d,u). From Eq. (fl7|) , we may read off the properties F(u,u) = —u 
and F(u,0) = C(2)u. For m D = 0, Eq. ([16]) coincides with Eq. (12) of Ref. 

Inserting Eqs. (|l~2]JT6|) into Eq. @, we obtain the general expression for the (U,D) con- 
tribution to T [H — > ff) at next-to- leading order, 
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For the reader's convenience, we list $teH[ for positive argument. For < r < 1, one has 
1 
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and, for r > 1, 
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where Cl 2 and CI3 are the (generalized) Clausen functions of second and third order [ p0|| . 
respectively, $ = arcsin-^/r, p± = y/r ± \/r — 1 , <p = lnp + , 7 = ln(p+ + p_), and x = 
ln(p + — p_). It is useful to know the expansions of H[ appropriate to the various limiting 
cases. They are 

13 1 22 

H[(r) = 6C(3) + 3C(2) - - + —r + 0(r 2 ) (21) 
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in the heavy-quark limit (r <C 1), 



H[(r) = C(2) (-12/1 - 6 In 2 + ^) - ^C(3) + ~ + 37r v / T^ : + - r)) (22) 

at threshold (r^l), and 

Ke^(r) = 6 7 2 - 3 7 - 6C(2) - ? + ^(-6 7 + 1) + ^\ (23) 

in the light-quark limit (r ^> 1). From the above results we can glean the leading QCD 
correction to Kf for f^b. Inserting Eq. into Eq. (|18D , with mu = m t ^> M#/2 and 
tub — 0, and comparing the result with Eq. ([I]), one finds the corrected value, 

K, = 7 - 3 (((2) + |) C F ^ 

- MM)* 

w 7- 4.004 a s , (24) 

where terms of (9 (M# / ) have been suppressed. We recover the notion that, in Electroweak 
Physics, the one- loop O (Gprn 2 ) terms get screened by their QCD corrections. 



3 Dependence on the quark-mass definition 

So far, we have employed the on-shell definition of quark mass, i.e., we have evaluated 
the counterterm diagrams of Fig. |l](b) and similar diagrams for the J^-boson self-energy 
using SrriQ in the form specified in Eq. ([11]) . This is certainly a reasonable choice. In the 
approximation of neglecting the p 2 dependence of the imaginary part of the quark self-energy, 
the on-shell mass coincides with the real part of the complex pole position, i.e., with the 
physical mass, which is a constant of nature p4| . In general, the physical mass is close 



to what is determined experimentally, e.g., in quarkonium spectroscopy. In the case of the 
top quark, it is approximately the physical mass that is being extracted at the Tevatron 
and will be at future e + e~ linear colliders, since, in the propagation of the t and i quarks 
between the production and decay vertices, configurations near the mass shell are greatly 
enhanced kinematically. As a matter of principle, however, this mass convention is arbitrary, 
and one might as well adopt another one. For, when the perturbation series is summed, the 
final result should not depend on the selected scheme. Yet, this holds no longer true when 
the perturbation series is truncated. In general, the finite-order results depend also on the 
renormalization scales of the quark masses. Scheme and typical scale variations may be used 
to estimate the theoretical uncertainty due to the unknown higher-order corrections. 

In perturbative-QCD calculations, the quark masses are frequently defined according to 
the MS scheme. In this case, 5rriQ collects just the pole in e, 

^ = -^C/-. (25) 
tuq Atx e 
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The relationship between the on-shell mass and the MS mass, rriQ, may be read off from 
Eq. (P Egg, 



m Q = m Q 
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The mass renormalization scale, n, must be chosen judiciously according to the problem at 
hand so as to minimize higher-order corrections. The O (a 2 ) term of Eq. ( |2"6Tj may be found 
in Ref. | 36fl , but we shall not need it here. 

In the following, we shall translate the results of Sect. 2 to the MS scheme. That is, 
we have to repeat the evaluation of the two-loop counterterm diagrams using SttTq instead 
of SrriQ. To this end, we may exploit the fact that the corresponding expressions may be 
constructed from the one- loop amplitudes by variation, e.g., 
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and similarly for the MS scheme. Consequently, the MS version of Tl HH may be obtained 
from Eq. (|12|) by including the term 
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Since 5mQ — 5mq is devoid of ultraviolet singularities, knowledge of the 0(e) term of n HH 
is not necessary. After carrying out this operation in Eqs. (§,0), we may represent the result 
by assigning shifts to the various items appearing in Eqs. (0,0), 
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(31) 

(U <-> D). (32) 



A complimentary set of shifts appropriate to the gauge-boson vacuum polarizations at arbi- 
trary four-momentum may be found in Ref. | 37| . Similarly to Eq. @, Eq. (^Ij) is valid for 
r < 0. To compute the shift of Eq. ([18]), we need A"ReH[ for r > 0. This is 



AH[(r) = (3£ Q - A) 

for < r < 1 and 

AUeH[{r) = {U Q - 4) 

for r > 1. 
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In the case of the (t, b) contribution to the fermionic decay widths of an intermediate- 
mass Higgs boson, we may set mjj = m t ^> M#/2 and mjj = 0. Then, the shift in Eq. (|T8| ) 
becomes 

4 2tt 2 V2 vr V ; V 6 V m * / / 

Thus, the MS version of Eq. fl2~3|) is given by 



^/ = 7+(-3C(2) + y-|^)^, (36) 

where terms of O (M H /mf) have been neglected. For \i = m t , this is 7 + (2/3)(19/7r — n)a s ~ 
7+ 1.938 a s . That is, the magnitude of the QCD correction is about half as large as in the 
on-shell case and its sign is opposite. 



4 Numerical analysis 



We are now in a position to explore the phenomenological consequences of our results. To 
start with, we specify the input values for our numerical analysis. We use Mw = 80.24 GeV 

1091 



91.1895 GeV 



4T|, and m t 



174 ± 16) GeV p 



polarization according to Ref. [|3 
at v/i = 91.175 GeV 
light quarks otherwise 



1.777 GeV |0|, m b = 4.72 GeV 
We parameterize the hadronic contribution to the photon vacuum 
with the updated reference value Aahadrons 



14 



0.0283 

and use an equivalent set of effective mass parameters for the 



(5) 

the asymptotic scale parameter appropriate to five active flavours, A^ 



For a s (/i), we employ the MS formula of Ref. j45| and fix 

by requiring that 

a s (Mz) = 0.124 |j39"| . Unless stated otherwise, we shall choose fi = m t for (t, b) contributions 
and /i = Mh else. All other input parameters are adopted from Ref. J|6] . 

In Fig. 0, we show versus Mh the radiative corrections to the leptonic decay widths of 
the Higgs boson originating from quark loops (dotted lines) and one-gluon exchanges within 
these loops for m t = (174 ± 16) GeV. These corrections are mainly due to the (t, b) doublet 
and do not depend on the flavour of the final-state fermions. As we have observed already 
in Sect. 2, the QCD correction reduces the one-loop quark contribution. For m t = 176 GeV, 
the screening effect ranges between 6.5% and 10.4% for Mh between 60 and 200 GeV. 

The full electroweak correction does depend on the produced flavour, as may be seen 
already from the approximation of Eq. ([!]). Figure presents the Mh dependence of the full 
electroweak one-loop |L8| (dotted lines) plus QCD two-loop correction to T(H — > t + t~) for 
m t = (174 ± 16) GeV. Since the QED correction depends logarithmically on Mh, the slopes 
of the curves are steeper than in Fig. 0. At one loop, there is a large cancellation between the 
bosonic and fermionic contributions |]I8fl . For m t = 174 GeV, their sum is in fact negative 
and relatively small, so that the QCD correction enhances significantly the size of the total 
correction, by 50.0% (9.7%) at M H = 60 GeV (150 GeV). 

Figures |2],|3| refer to the on-shell definition of quark mass. In Figs. |]|| we study how 
the radiative corrections to the leptonic decay widths of the Higgs boson are affected when 
the quark mass renormalization is converted from the on-shell scheme to the MS scheme. 
For a meaningful study of the scale dependence, it is necessary to distinguish between the 
renormalization scales of the quark mass and the strong coupling constant, fi m and /i c , 
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respectively. The scale /i that occurs explicitly in the above formulae must be identified 
with ji m . Since the quark contributions scale like m 2 [see Eqs. (0,0)], we may restrict 
our considerations to the top quark and neglect the masses of the other quarks. Figure [| 
visualizes the one- and two-loop evaluations in the on-shell and MS schemes as a function of 
the top-quark on-shell mass, mt, for Mh = 100 GeV. The solid and dashed lines represent 
the one- and two-loop on-shell results, respectively. The dot-dashed line corresponds to the 
one- loop evaluation with m t replaced by m t as given by Eq. (|26|) with fi m = m t . Inclusion 
of 5i + A5\, evaluated with fn t , then leads to the dotted line. The difference between the 
dashed and the dotted lines is hardly visible on the plot, which indicates that the scheme 
dependence at two loops and beyond is negligibly small. 

In addition to the trivial \x c dependence via a s , which is present also in the on-shell 
scheme, the MS analysis depends on the mass renormalization scale, fi m , too. A plausible 
choice is /i m = m t , since this eliminates the terms proportional to if By the same token, it 
seems unnatural to choose \x m very different from m t , since this renders £ t artificially large. 
In Fig. H, we analyze the scale dependence of the two-loop MS result for Mh = 100 GeV 
and m t = 174 GeV. As usual, we introduce a dimensionless scale parameter, £, which we 
vary from 1/4 to 4. The dotted line corresponds to the choice fi c = £m 4 , \x m = m t , which 
could be studied also in the on-shell scheme. Since the QCD corrections are known only 
to leading order, the £ dependence is monotonic. It is expected to flatten when three-loop 
QCD corrections will be taken into account. The dashed line represents the evaluation with 
fi c = m t , \x m — C m t- Here, we observe a stabilization already at two loops. This is due to a 
cancellation between the \i m dependence induced in the one-loop result via m t and the one 
carried by the genuine two-loop term. There is no obvious reason to choose /i m different 
from /i c . When we identify the two scales and set /i c = /i m = £m f , we obtain the solid 
curve, which assumes its maximum very close to one, at £ = 0.958. This choice may thus be 
advocated by appealing to the principle of minimal sensitivity HTI . 



5 Conclusions 

In the Standard Model, low- and intermediate-mass Higgs bosons decay preferably into 
fermion pairs, and the one-loop electroweak corrections to the respective decay rates are 
dominated by terms of 0{Gpm 2 ). The leading two-loop corrections to these terms are of 
O (asGprnl) and O {G 2 F m t ). We have calculated the two-loop gluon correction to the shift 
in the fermionic decay rates of the Higgs boson induced by pairs of virtual quarks with 
arbitrary masses. This correction is of O (a s GVm;Q, where mg is the mass of the heaviest 
quark, and does not depend on the produced fermion flavour. In the case of lepton-pair 
production, this is the only QCD correction up to two loops. In the case of the hadronic 
decays of the Higgs boson, this is the only source of O (a s Gpm 2 ) corrections. A special 
situation arises in the case of bb final states. Here, the universal O {Gf^i 2 ) term of the 
electroweak one-loop correction is almost entirely cancelled by similar contributions to the 
bbH vertex and the 6-quark wave function induced by the charged current. The two- loop 
O (a s GFm 2 ) corrections to the latter contributions are currently under study pE8"[ . 

For M H = 60 GeV (200 GeV), QCD effects screen the positive universal O (Gprn 2 ) term 
by 6.5% (10.4%), provided the top quark mass is renormalized according to the on-shell 
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scheme. In the MS scheme, the QCD correction has a different sign and its magnitude 
is roughly half as large. However, the sum of the one- and two-loop corrections is very 
insensitive to the choice of mass renormalization scheme. In the MS evaluation with a 
common renormalization scale, /i, for the strong coupling constant and the top-quark mass, 
the total correction is stable under variations of \x in the vicinity of \i = m%. This indicates 
that the effect of QCD corrections beyond two loops is likely to be insignificant. 
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FIGURE CAPTIONS 



Figure 1: Feynman diagrams pertinent to the QCD corrections to the Higgs-boson self- 
energy: (a) gluon exchanges and (b) counterterms. 

Figure 2: 0(a) and 0(aa s ) radiative corrections to T (H — > £ + £~), with £ = e,[i,r, due to 
virtual quarks as a function of M H for m t = (174 ± 16) GeV. 

Figure 3: Full 0(a) and 0(aa s ) radiative corrections to T (H — > t + t~) as a function of M# 
for m t = (174 ± 16) GeV. 

Figure 4: 0(a) and 0(aa s ) radiative corrections to T (if — > with £ = e,fi,T, due to 

virtual quarks as a function of m t for M# = 100 GeV. The evaluations using the on-shell 
and MS definitions of the top-quark mass are compared. All other quark masses are set to 
zero. 

Figure 5: 0(aa s ) radiative corrections to T (H — > £ + £~), with £ — e, /i, r, for Mh = 100 GeV 
and m t = 174 GeV evaluated in the MS scheme with \i c = £,m t , fi m = m t (dotted line), 
fi c = m t , /i m = £m t (dashed line), and /i c = ji m = ^m t (solid line). All other quark masses 
are set to zero. 



14 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405299vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405299vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405299vl 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405299vl 



